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Abstract. After 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March 
that lasted several months. During the active phase, several bursts were detected, and many aspects of the X-ray emission 
changed. We report on the discovery of six X-ray bursts, the first ever seen from this AXP in ~10 years of Rossi X-ray 
Timing Explorer monitoring. All the bursts occurred in the interval between 2006 April 6 and 2007 February 7. The burst 
durations ranged from 8— 3x 10 3 s as characterized by Tgo. These are very long durations even when compared to the broad 
Tqo distributions of other bursts from AXPs and Soft Gamma Repeaters (SGRs). The first five burst spectra are well modeled by 
simple blackbodies, with temperature kT ~ 2 — 6 keV. However, the sixth and most energetic burst had a complicated spectrum 
consisting of at least three emission lines with possible additional emission and absorption lines. The most significant feature 
was at ~ 14 keV. Similar 14-keV spectral features were seen in bursts from AXPs IE 1048.1-5937 and XTE J1810-197. If 
this feature is interpreted as a proton cyclotron line, then it supports the existence of a magnetar-strength field for these AXPs. 
Several of the bursts were accompanied by a short-term pulsed flux enhancement. We discuss these events in the context of 
the magnetar model. 
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INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) are isolated neutron 
stars that show pulsations in the narrow range of 2- 
12 s. Their observed 2-10 keV X-ray luminosities (~ 
10 33 — 10 35 erg s _1 ) cannot be accounted for by their 
available spin-down energy. It is widely accepted that 
AXPs are magnetars - young isolated neutron stars pow- 
ered by their high magnetic fields The inferred 
surface dipolar magnetic fields of AXPs are all above 
5.9 xlO 13 G. The magnetar model was first proposed to 
explain Soft Gamma Repeaters (SGRs). SGRs show per- 
sistent properties similar to AXPs, but they were first 
discovered by their enormous bursts of soft gamma rays 
(> 10 44 erg) and their much more frequent, shorter, and 
thus less energetic bursts of hard X-rays. To date, SGR- 
like X-ray bursts have been observed from four AXPs, 
thus solidifying the connection between the two source 
classes |Hl|3,|5[0]. For a review of magnetar candidates 
see Woods and Thompson 10] . 

Thus far, only the magnetar model can explain the 
bursts observed from SGRs and AXPs [1]. The internal 
magnetic field exerts stresses on the crust which can lead 
to large scale rearrangements of the external field, which 
we observe as giant flares. If the stress is more localized, 
then it can fracture the crust and displace the footpoints 
of the external magnetic field which results in short X- 
ray bursts. The highly twisted internal magnetic field also 



slowly twists up the external field and it is believed that 
the magnetospheres of magnetars are globally twisted 
01 . Reconnection in this globally twisted magnetosphere 
has also been proposed as an additional mechanism for 
the short bursts flSj] • 

In addition to bursts, AXPs and SGRs exhibit pulsed 
and persistent flux variations on several timescales. An 
hours-long increase in the pulsed flux has been seen to 
follow a burst in AXP IE 1048.1-5937 (H. On longer 
timescales, AXPs can exhibit abrupt increases in flux 
which decay on ^week-month timescales. These occur 
in conjunction with bursts and are thought to be due to 
thermal radiation from the stellar surface after the depo- 
sition of heat from bursts. Such flux enhancements have 
been observed in SGRs (see Woods et al. Ill ill for exam- 
ple). The flux enhancement of AXP IE 2259+586 dur- 
ing its 2002 outburst was also interpreted as burst after- 
glow lfl2ll . however, a magnetospheric interpretation has 
also been proposed fH]. AXP IE 1048.1-5937 exhibited 
three unusual flux flares. In the first two, the pulsed flux 
rose on week-long timescales and subsequently decayed 
back on time scales of months lfl4l[l5ll . Although small 
bursts sometimes occur during these events [10], burst af- 
terglow cannot explain the flaring, thus these variations 
have been attributed to twists implanted in the external 
magnetosphere from stresses on the crust imposed by the 
internal magnetic field. AXPs XTE J1810-197 and the 
AXP candidate AX J 1845-0258 have also exhibited flux 



variations, however it is not clear whether these were of 
the abrupt rise type as in IE 2259+586 or the slow-rise 
type as in IE 1048.1-5937. Finally, AXP 4U 0142+61 
has exhibited the longest timescale flux variations, in 
which the pulsed flux increased by 19 ±9% over a period 
of 2.6 years fH]. 



ANALYSIS AND RESULTS 

All data presented here are from the Proportional 
Counter Array (PCA) aboard the Rossi X-ray Timing 
Explorer (RXTE). The PCA is made up of five identical 
and independent proportional counter units (PCUs). 
Each PCU is a Xenon/methane proportional counter 
with a propane veto layer. The data were collected in 
either GoodXenonwithPropane or GoodXenon 
mode which record photon arrival times with — 1-jiis 
resolution and bins them with 256 spectral channels in 
the -2-60 keV band. 



Burst Analysis 

We have been monitoring 4U 0142+61 with the PCA 
for nearly a decade. Currently, it is observed bi-monthly 
with a typical observation length of 5 ks. For each 
monitoring observation of 4U 0142+61, using software 
that can handle the raw telemetry data, we generated 
31.25 ms lightcurves using all Xenon layers and only 
events in the 2-20 keV band. These lightcurves were 
searched for bursts using our automated burst search al- 
gorithm introduced in Gavriil et al. [3] and discussed fur- 
ther in Gavriil et al. [17]. In an observation on 2006 April 
6, we detected a significant burst, and four more bursts 
were detected in a single observation on 2006 June 25. 
The sixth and most energetic burst was detected on 2007 
February 7. There were 3, 3, and 2 PCUs on at the times 
of the bursts for the April, June and February observa- 
tions, respectively. The bursts were significant in each 
active PCU. 

To further analyze these bursts we created event lists 
in FITSQ format using the standard FTOOLS0 For con- 
sistency with previous analysis of SGR/AXP bursts we 
extracted events in the 2-60 ke V band. These events were 
barycentered using the position found by Patel et al. Ill 811 
for the source. The burst lightcurves are displayed in 

Fig.m 

Before measuring any burst parameters we deter- 
mined the instrumental background using the FTOOL 
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FIGURE 1. The histograms are the 2-60 keV burst 
lightcurves binned with 1/32 s resolution as observed by RXTE. 
The thick curves are the best fit exponential rise and exponen- 
tial decay model. 



pcabackest. We extracted a background model 
lightcurve using the appropriate energy band and 
number of PCUs. pcabackest only determines the 
background on 16 s time intervals, so we interpolated 
these values by fitting a polynomial of order 6 to the 
entire observation, which yielded a good fit for each 
observation. 

The burst peak time, rise time and peak flux were de- 
termined using the methods described in Gavriil et al. 
ifr/Tl . Usually, to measure the fluence for SGR and AXP 
bursts, we subtract the instrumental background for the 
lightcurve, integrate the light curve and fit it to a step 
function with a linear term whose slope is the "local" 
background rate. The fluence in this case is the height of 
the step function. Although this technique worked well 
for the first burst, which was a short isolated event, it 
was not appropriate for bursts 2, 3, and 4 because they 
had overlapping tails, and bursts 5 and 6 had tails that 
extended beyond the end of the observation. Thus, we 
opted to fit the bursts to exponential rises with decay- 
ing tails. Our model fits are overplotted on the bursts in 
Fig. Q] As is done for y-ray bursts and SGR and AXP 
bursts, we characterized the burst duration by Tyo, the 
time from when 5% to 95% of the total burst counts have 
been collected. To determine the T90 duration we inte- 
grated our burst model and numerically determined the 
5% and 95% time crossings. All burst temporal parame- 
ters are presented in TableQ] 

Burst spectra were extracted using all the counts 
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within their 790 interval. Background intervals were ex- 
tracted from long, hand-selected intervals prior to the 
bursts. Response matrices were created using the FT00L 
pcarsp. The burst spectra were grouped such that there 
were at least 15 counts per bin after background sub- 
traction. Burst spectra, background spectra, and response 
matrices were then read into the spectral fitting pack- 
age XSPEC0 vl2.3.1. The spectra were fit to photoelec- 
trically absorbed blackbodies using the column density 
found by Durant and van Kerkwijk lfl9ll . Only bins in 
the 2-30 keV band were included in the fits. The black- 
body model provided an adequate fit for bursts 1 through 
5. burst 6, however, was not well modeled by any sim- 
ple continuum model because of the presence of emis- 
sion lines (see Fig|2]panels 1A and IB). These features 
showed clear temporal variability but they were most 
prominent near the onset of the burst (see Fig. [3]). 



Pulsed Flux Analysis 

For each of the three observations containing bursts, 
we made two barycentered time series in count rate 
per PCU, one for the 2—4 keV band and the other for 
4—20 keV. We only included the photons detected by 
the PCUs that were on for the entire duration of the 
observation. The time resolution was 1/32 s. We removed 
the 4 s centered on each burst from each time series. 
Then, we broke each time series into segments of length 
^500 s. For each segment, we calculated the pulsed flux 
using two different methods. 

First, we calculated the RMS pulsed flux using the 
Fourier decomposition method described by Woods et al. 
UH], only incorporating the contribution of the first 5 
harmonics for consistency with lfl6ll and j20Tl . While 
least sensitive to noise, the RMS method returns a pulsed 
flux number that is affected by pulse profile changes 
(Archibald et al. in prep.). So to confirm our pulsed flux 
results, we also used an area-based estimator to calculate 

the pulsed flux, PF = «o ^W^' wnere a o — wllf=iPh 

i refers to the phase bin, N is the total number of phase 
bins, pi is the count rate in the i" 1 phase bin of the pulse 
profile, and p m [ n is the average count rate in the off -pulse 
phase bins of the profile, determined by cross-correlating 
with a high signal-to-noise template, and calculated in 
the Fourier domain after truncating the Fourier series to 
5 harmonics. The results are shown in Figure|4] Note the 
significant increase in the 4—20 keV pulsed flux in the 
2007 June observation following the cluster of bursts. 
This increase is not present in 2—4 keV. Also note the 
significant rise and subsequent decay of the pulsed flux 
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FIGURE 2. Burst spectra of all AXP bursts with significant 
emission lines as observed by RXTE. (1A) Burst spectrum of 
4U 0142+61 burst 6. The dotted line indicates the continuum 
(blackbody) component of the best fit model. IB: Residuals af- 
ter subtracting the continuum component of the best fit model. 
2A and 2B: Same but for burst 4 of XTE J1810-197 [see@]. 
3A and 3B: Same, but for burst 3 of IE 1048.1-5937 [seefp]. 
4A and 4B: Same, but for burst 1 of IE 1048.1-5937 [seefaT 



following the large 2007 February burst. 



DISCUSSION 

We have discovered six bursts from AXP 4U 0142+61. 
These bursts all occurred between 2006 April and 2007 
February, and were the only ones ever observed from 




FIGURE 3. Dynamic spectrum of burst 6. The wedge indi- 
cates the number of counts as a function of time and energy. 
Notice how at later times the contribution of the features at 
~4 keV and ~8 keV exceed that of the 14 keV feature. 
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FIGURE 4. RMS and area pulsed flux within the observa- 
tions containing bursts. Each column corresponds to one ob- 
servation. In each column we have, descending vertically, the 
1-s resolution lightcurve with the bursts indicated, the 2^1 keV 
RMS pulsed flux, the 2^1 keV area pulsed flux, the 4-20 keV 
RMS pulsed flux, and the 4-20 keV area pulsed flux. The dotted 
line in each of the pulsed flux plots shows the average of the 
pulsed fluxes obtained after segmenting and analyzing the time 
series of the observation immediately prior to the one shown. 



this source in ~10 years of monitoring. After the first 
burst 4U 0142+61 exhibited a timing anomaly and pulse 
profile variations (Gavriil, Dib, & Kaspi in preparation). 
Together with the short-term pulsed flux increase, the 
simultaneity of all these phenomena clearly identifies 
4U 0142+61 as the origin of the bursts. 

Woods et al. JH first argued that there appear to be 
two classes of magnetar bursts. Type A bursts are short, 
symmetric, and occur uniformly in pulse phase. Type B 
bursts have long tails, thermal spectra, and occur prefer- 



entially at pulse maximum. Woods et al. ^ noted that 
type A bursts occur predominately in SGRs and type B 
bursts occur predominately in AXPs, and this was af- 
firmed by Gavriil et al. II 1 CXI . Both argue that Type A and 
Type B bursts are produced by different mechanisms. In 
the magnetar model bursts can either be due to the rear- 
rangement of magnetic field lines anchored to the surface 
after a crustal fracture 1 1], or due to reconnection in the 
upper magnetosphere [9]. Woods et al. 15J and Gavriil 
et al. [10] argue that Type B bursts are due to the for- 
mer and Type A bursts are due to the latter. The bursts 
reported here all had very long tails, 790 > 8 s, suggest- 
ing they are of Type B. However, the bursts did not occur 
preferentially at pulse maxima. 

The line-rich spectrum of burst 6 is intriguing. Three 
significant features are seen at ^4, ^8 and ^14 keV. 
The most significant emission feature at ^14 keV is par- 
ticularly interesting. Emission features at similar ener- 
gies were observed from two out of the three bursts from 
IE 1048. 1-5937 |f3l [Toh and in one out of the four bursts 
from XTE 11810-1970, H. We have reanalyzed these 
burst spectra in a consistent manner as for 4U 0142+61. 
In Fig.|2]we plot the spectra of all AXP bursts with emis- 
sion lines in their spectra. Notice that all the spectra have 
features that occur between 13 and 14 keV and are very 
broad. There is clear evidence for features at ~4 and 
~8 keV in the 4U 0142+61 burst; however, note that 
there is subtle evidence for these features in some of the 
other burst spectra as well. 

If the ~+4 keV feature is interpreted as a proton cy- 
clotron feature then we can infer the surface magnetic 
field strength of the star. For a line of energy E the mag- 
netic field strength is given by 

Setting m equal to the proton mass we obtain B — 2.2 x 
10 15 (E/14- keV) G. This field estimate is much greater 
than that derived from the spin down of the source, how- 
ever the burst spectroscopic method measures the field at 
the surface which can be multipolar, while the spin down 
measurement is sensitive to the dipolar component. 

The feature also could be an electron cyclotron fea- 
ture at the surface. Replacing m with the electron mass in 
Eq.Q]we obtain B= 1.2 x 10 12 (£/14 keV) G. This field 
is two orders of magnitude less than the spin down field. 
However, if the feature occurred higher up in the magne- 
tosphere then the field would be greatly reduced. Thus, 
an electron cyclotron feature from a burst which occurred 
in the upper magnetosphere cannot be precluded. 

Although these features can, in principle, be pro- 
ton/electron cyclotron features there are many problems 
with interpreting them as such. First, it is not clear 
why three different sources, with different magnetic field 
strengths would exhibit features with similar energies. 



Moreover, it is unclear why these features are seldom 
seen and have not been seen in other high signal-to- 
noise bursts. Detailed modeling of burst spectra is def- 
initely warranted. The fact that these features occur at 
similar energies, despite having different magnetic field 
strengths, may suggest that they are actually atomic lines, 
and could possibly provide new insights into the com- 
position of the crust and atmosphere of these enigmatic 
objects. 
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